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Abstract 

The survey reports theoretical studies involving positronium (Ps) - atom scat- 
tering. Investigations carried out in last few decades have been briefly reviewed in 
this article. A brief description of close-coupling approximation (CCA), the first- 
Born approximation (FBA) and the Born-Oppenheimer approximation (BOA) for 
Ps-Atom systems are made. The CCA codes of Ray et al [1-6] are reinvestigated 
using very fine mesh-points to search for resonances. The article advocates the 
need for an extended basis set & a systematic study using CCAs. 


1 Introduction 

The discovery of positronium (Ps), a hydrogen (H) like atom formed by a positron 
(e + ) and an electron (e~) in 1953 by Deutsch [7], is an invaluable achievement 
of the modern science. The atom is itself its anti-atom. Interesting property is 
that its charge and mass centers coincide. The existence of Ps was first predicted 
by S. Mohorovicic in 1934 [8] and later on by Ruark in 1945 [9] based on their 
theoretical investigations which encouraged Deutsch in 1949 [10] to carry out his 
experimental investigation. The concept of anti-particle was introduced by Dirac 
in 1928 [11-12] as an anti-electron which was later experimentally observed by 
Anderson in 1932 [13] and was named as positron. 

The collision physics is the most important area of the modern science. The 
most successful atomic-model by Bohr (1913) and Sommerfeld (1916) was based 
on the well-known Rutherford scattering (1911) experiment and theory; needless 
to say that the atomic concept was conceived more than a century ago by Dalton 
(1808), Gay-Lussac (1808), Avogadro (1811), Maxwell (1811), Mendeleef (1869) 
and others. 

The theoretical studies on Ps and gas atom scattering was initiated by Massey 
and Mohr [14-15] in 1954 due to its interesting properties. The present progress 
is due to the continued interests of Fraser [16-18], Fraser & Kraidy [19], Hara 
& Fraser [20], Martin h Fraser [21], Barker & Bransden [22-23], Bransden [24], 
Drachman & Houston [25-26], Drachman [27-29], Au & Drachman [30], DiRienzi 
h Drachman [31-32], Schrader et al [33] and the realization of the fact of all 
the recent workers doing theories [1-6,34-51] and experiments [52-57] in different 
countries and different places of the world. The recent news published on ’Universe 
Today’ dated June 8, 2005 has the headline, ’’The search for Positronium”. This 
could be an important source of many new ideas leading to new physics. 

The positrons and Ps are currently employed in the exploration of fundamental 
effects ranging from condensed matter physics to astrophysics as well as in the di- 
agnostics of living biological systems and of the electronic and structural properties 
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of industrially important materials. The study of collisions of neutral positronium 
with other atoms (and molecules) has a number of novel features arising from its 
light mass and from the coincidence of its centres of charge and mass. Aside from 
this intrinsic interest, Ps collision data are expected to enhance our understanding 
of positron slowing down in dense media, Ps diffraction from surfaces, interactions 
of Ps injected into a plasma, etc. Various industrial applications of positron and 
positronium beams are discussed by Coleman [58]. 

Two interacting atoms mutually induce a symmetrical pair of dipoles [59]. 
The induced polarization potential is the same for electron-atom and positron- 
atom scattering, in both cases they are attractive [24]. Ps is a highly polarizable 
atom, its polarizability is eight times higher than that of H. They can exist in two 
different spin states depending on the spin of positron and electron. The singlet 
is known as para-Ps with a life time ~ 10 -10 s and the triplet is known as ortho- 
Ps having life time ~ 10 -7 s. The electrons are indistinguishable particles. Two 
electrons can interchange their positions, the phenomenon is known as exchange. 
The exchange is highly important at low incident energies in presence of more than 
one electron. The electron spins can be up or down. If both the spins are parallel 
and total spin is 1, it is known as triplet (-) state; if antiparallel, the total spin is 
0 and is known as singlet (+) state. 

2 Resonances 

In collision physics, the existence of a resonance is an important phenomenon. 
When a microscopic moving object which is a wave, enters into the scattering 
chamber near the target, it faces interactions. When it comes out of the scattering 
zone, the original incident wave gathers a phase shift and the new wave is known 
as scattered wave. The change in phase which is named as phase shift is the 
parameter that carries the information of the scattering process. A rapid change 
in phase shift by 7r radian in a very narrow energy interval of the incident wave is 
an indication of the presence of a resonance. It indicates the existence of a bound 
system if in the s-wave elastic scattering and below threshold of excitation. One 
can calculate the width of a resonance to get the life time (r) of the newly formed 
system. 

It is an extremely difficult job to detect a resonance since successful identifi- 
cation needs (i) a very accurate calculation and (ii) sufficient computation facil- 
ities. A large number of mesh-points in a very small energy interval, generally 
~ 10 -2 — 10~ 3 eV is required. It necessitates a high-speed computer with a suf- 
ficient memory, the knowledge of mathematical computation and programming 
languages e.g. FORTRAN. 

We are interested to discuss the scattering processes at low energies. At high 
energies the projectile usually ignores all the important delicate interactions with 
target, so generally it carries no valuable information. At very low energies below 
excitation threshold, elastic scattering is the only real process. But from just above 
the threshold different excitation, ionization etc. channels start to open. It is a 
very difficult task to study the scattering processes at intermediate energies due to 
presence of many different channels and a very close-coupling among them. The 
total cross section is the sum of the integrated cross sections of all these channels. 
Besides these, the partial wave contribution from higher angular momenta start 
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to dominate above the threshold. Many close peaks may arise in the cross section 
at the intermediate energies due to opening of different channels and their mixing. 
These are also named as resonances and have been studied by Higgins & Burke 
[60-61], Sarkar et al [62-63]. However Zhou et al [64] commented against such 
resonances. So the subject of above threshold resonances needs more investigation. 
Above threshold, the phase shift becomes a complex quantity and is known as eigen 
phase shift. 


2.1 Below threshold 

The kind of resonances we are interested in, is completely different from the above 
threshold resonances. This type of a resonance is feasible only below threshold of 
first excitation when elastic cross section is the total cross section and only s-wave 
dominates. One needs to investigate both the phase shift and the cross section, but 
should be careful to use the formulation which derives the cross section directly 
from amplitude and not from phase shift. If it is a true resonance, this should 
be reflected both in phase shift and in cross section. It is actually to take a 
precautionery measure for successful detection since the calculation of phase shift 
involves a tan-inverse function which may create a numerical error. 

The phase shift Si can be decomposed as 
Si = + m- 

£i corresponds to the hard sphere scattering or non-resonant part; it does not 
depend on the shape and depth of the potential. The term rji depends on the 
details of the potential. The quantities £/ and r]i vary, in general, slowly and 
smoothly with the incident particle energy. But in certain cases r]i may vary 
rapidly in a small energy interval of width T about a given energy value Er such 
that we can write 

»,=„« = ton ' 1 2(£r _ e) 


In that energy interval the phase shift is therefore given approximately by 

h - & + »?/*• 

The physical significance of a narrow resonance can be inferred by examining 
the amplitude of the radial wave function inside the interaction region. The prob- 
ability of finding the scattered particle within the potential is much higher near 
the resonance energy E = Er, so that in that case the particle is nearly bound 
in the well. Thus the resonance may be considered as a metastable state whose 
lifetime r, which is much longer than a typical collision time, can be related to 
the resonance width T by using the uncertainty relation AtAE > h. Thus, with 
A t ~ t and A E ~ T, we have r ~ 

The shape of the cross section curve near a resonance as a function of energy 
depends on the non-resonant phase shift £/. For the s-wave scattering it is 

_ sin 2 ti(E R - E) 2 + cos 2 ZiT 2 + sin2ji{E R - E)T/2 
° l ~ (Er -E)*+¥- ' ~ 

Two limiting cases for non-resonant phase shift are 0 and ir/2. In the first case 
the above equation becomes 
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r 2 

° l ~ {E R -EY + V i 

which is symmetric and represents a rise in cross section at the resonance energy. 
In the other case 


(Er - E ) 2 

^ (Er - Ef + £ 

which is also symmetric but goes down to zero at the resonance energy. If the 
non-resonant phase shift gets some other value then all sorts of forms of the cross 
section can occur. 

Resonance in the singlet channel in positronium (Ps) and hydrogen (H) scatter- 
ing was reported by many workers using different approaches [25,35,37-38,65-67]. 
It was first predicted by Drachman et al [25] using a Feshbach formalism with sta- 
bilization and complex rotation methods. Such a resonance using CCA schemes 
is reported for the first time by Ray [1-2,68-72]. 

3 Theory 
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Figure 1: The s-wave elastic phase shifts below inelastic thresholds in Ps-H scattering 
using target-inelastic CCA theory. 


3.1 Close-coupling approximation (CCA) theory 

The close-coupling approximation (CCA) is a successful theory to study the low 
energy scattering phenomenon in atomic physics. The formulation was given by 
Massey [73]; he is known as the father of atomic physics. Massey applied the 
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Figure 2: The s-wave elastic cross sections below inelastic thresholds in Ps-H scattering 
using target-inelastic CCA theory. 

theory in e~ - atom scattering. Later Burke et al [74] successfully used this theory 
for e + - atom scattering. Now a days many different groups in the world are using 
this theory for e~/e + - atom scattering. It was Fraser [16] who used it first for 
Ps-H scattering. Ray & Ghosh [47-48] merits the credit because they are the 
first who supplied detailed and converged results. They used a momentum space 
formalism introduced by Calcutta group [75] whereas Fraser used a coordinate 
space formalism to write the coupled integral equations. They again add more 
channels in the CCA basis [1-6,68-72,76-78]. The studies of Fraser [16-18], Fraser 
et al [19], Hara et al [20] were confined to static-exchange model i.e. considering 
only the elastic channel in the basis and the H and He targets. Ray [4-6] extended 
the CCA theory in Ps and lithium (Li) scattering using the static-exchange and a 
two-channel CCA models. 

The theory is based on the very basic principle of quantum mechanics i.e. the 
eigen state expansion (ESE) methodology in which the total wave function of a 
quantum mechanical system is expressed as a linear combination of all possible 
states known as basis set. So one has to use a wide channel space, but practically 
it is not possible. The number of unknowns exceed the number of equations when 
non-spherical orbitals like p- and d- states of an atom are considered in channel 
space. So arose the necessity of an approximation. We should conserve the total 
angular momentum quantum numbers e.g. and ‘M‘. It makes the equations 
closed i.e. the number of unknowns are equal to the number of equations. This is 
known as CCA. The accuracy of the method depends on the choice of basis set. 

The total wavefunction of Ps-H system 4/^ satisfying the Schrodinger equation: 

^' ± (rp,ri,r 2 ) = ■E® ± (r p ,ri,r 2 ) (1) 
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Figure 3: The s-wave elastic phase shifts below inelastic thresholds in Ps-H scattering 
using projectile-inelastic CCA theory. 


is expressed as 

^ ± ( r p,n,r 2 ) = -^=(1±P 12 ) 
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Here P 12 stands for the exchange operator and Ri = ^(r p + rj) and pi — r p — r;; 
i=l,2; ri and r 2 are the position vectors of the electrons belonging to Ps and H 
respectively and r p , is that of the positron with respect to the center of mass of the 
system. U nt i t (r)/r and V n P i p (p)/p are the radial parts of the wavefunctions of H 
and Ps respectively and Fr 0 r(k, k', R)/R is the radial part of the continuum wave 
function of the moving Ps atom; To indicates all the quantities ntltn v l v LJ\JM of 
T at the initial channel. 

Projecting the Schrodinger Eqn.(l) just like the Hartree-Fock variational ap- 
proach and integrating over the desired coordinates, we can get a set of integro- 
differential equations which can be transformed into the integral equations like 
Lippmann- Schwinger applying the asymptotic boundary conditions. These cou- 
pled integral equations can be formed either in momentum space or in configuration 
space. We have used momentum space formalism [75]. The set of coupled integral 
equations obtained for the scattering amplitudes are as follows: 

± 


fnr, nl(k'.k) = + A Ej A " 


Dill 


( k/ » k/ 0/n''f'>i( k "’ k ) 


k" 2 - k 2 ,,^ + ie 
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Figure 4: The s-wave elastic cross sections below inelastic thresholds in Ps-H scattering 
using projectile-inelastic CCA theory. 

B ± indicate the Born-Oppenheimer [79] scattering amplitudes, plus(-h) is for 
the singlet channel and minus(-) for triplet channel. The formulation for the FBA 
and BOA matrix elements are discussed [80-86] and briefly described here. Simi- 
larly indicate the unknown scattering amplitudes for the singlet and the triplet 
channels respectively. The summation over n n l n is to include various channels. 

These three dimensional coupled integral equations involving f n 'i', n i{ k',k) and 
^n , z , ,nz(k , ,k) can be reduced to the corresponding one dimensional forms through 
partial wave analysis using the expansion like: 

= E E E E ( L 'l' r M' L m! r \J[M[) 

J,M L,Ml L' ,M' l Ji,Mi 


(Ll p , M L m p \J\M\)(j[l' t , M[m\\JM )(J 1 l t , M imt \JM) 


Yvm v (*) Ylml (k)r J± (k'n' p n' t l' p l' t , kn p n t l p l t ) (5) 

The resulting one dimensional coupled integral equations can be written as a 
matrix equation like 


[A]nxN Mnxi — PInxi (6) 

which can be solved by matrix inversion method. Here [A] is the scattering matrix 
of NxN type formed by Born and Born-Oppenheimer scattering amplitudes, [X] 
is the column matrix of N x 1 type formed by unknown CCA scattering amplitudes 
and [B] is again a column matrix of N x 1 type formed by the Born and Born- 
Oppenheimer [79] scattering amplitudes; the dimension(N) depends on the number 
of channels included in the expansion basis. We have calculated all the Born and 
Born-Oppenheimer amplitudes exactly following an analytic approach and then 
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Figure 5: The s-wave elastic phase shifts below inelastic thresholds in Ps-Li scattering 
using target-inelastic CCA theory. 

the remaining part is carried out following numerical approach using computer 
and FORTRAN programming. The two sets of one dimensional coupled integral 
equations of scattering amplitudes in momentum space for the singlet(+) and 
triplet (-) channels respectively, are solved separately for each partial wave(L). 

We employ different projectile-elastic and projectile- inelastic CCA-schemes to 
investigate Ps-H and Ps-Li scattering in the energy region below inelastic thresh- 
old. We study s-wave elastic phase shifts and s-wave elastic cross sections using 
a very fine mesh-points. We perform exact calculations for all the direct and ex- 
change matrix elements considering all the possible Coulomb interactions where 
the direct first-Born amplitudes vanish if the parity of Ps remains unaltered. 


3.2 First-Born approximation (FBA) 

The FBA amplitude for the scattering of the Ps by atomic target is expressed as, 
F f f = ~£j e- ikf - (x+ri/2 ) ^(|x-r 1 |)0}(r 2 ,r 3 ,..,r N )[F int ] 

<M r 2 ,r 3 , ..,r N )77i(|x - ri|)e lki (x+ri/2) dxdridr 2 <ir 3 ..dr N (7) 

where fj> is the reduced mass of the system; x, ri and rj are the position vectors 
for the positron, the electron in the Ps and the i-th electron in the target atomic 
system respectively w.r.t. center of mass of the system. N is the number of active 
electrons present in the system, rji and 77 / are respectively the initial and final 
state wavefunctions for the positronium atom whereas fa and (f>f are the initial 
and final state wavefunctions of the atomic target. V{ n t is the interaction potential 
due to electrostatic Coulomb interaction between the two atomic systems and is 
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Figure 6: The s-wave elastic cross sections below inelastic thresholds in Ps-Li scattering 
using target-inelastic CCA theory. 


given by 


_ Z Z ^ 1 ^ 1_ 

mt ~ x r\ |x - ri| + |n - rj 


( 8 ) 


i=2 1 ‘*1 i=2 1 

Making the substitution p = x — ri and R = (x + ri )/2 and performing the 
integrations, the above scattering amplitudes reduce to a general form like : 


Ffi — — {I Ps } {I target} 


(9) 


Here q = kj — kf , is the momentum transfer; kj and kf are the momenta of the 
projectile in the initial and final states respectively. Ip s and Itarget ? the form 
factors for Ps and target respectively are given by 


Ip, = I n}(p){e- i<i P/2 - e*-f>t 2 }rH{p)dp 


and 


(10) 


Itarget = [ <£/(r 2 , r 3 , r N )[Z - ^ e tq - ri> % (r 2 , r 3 , r N ) jj dr N ( 11 ) 
J i—2 i = 2 

Here Z represent the nuclear charge of target atom. 

3.3 Born-Oppenheimer approximation (BOA) 

The differential cross section in Ps-H scattering is defined as 


*_£ { j„ +0 p + i„- OP} 


(12) 
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if F and G are respectively the direct and exchange matrix elements in space part 
only. 

Accordingly the integrated or the total cross section (a) can be defined as: 


a(Ei) = ^-J dk f [\ | F(k f ) + G(k f ) | 2 + 1 | F(k f ) - G(k f ) | 2 ] (13) 

Where 

F (kf) = — J e- lkt - ni T]}(p 1 )^}{r 2 }[V F t }e iki - ni r] i (p 1 )^i{r2}dy:dr 1 dr2 (14) 

G{kf) = ~J e_lkf R2 ^/(^2)^/{ ri }[ Vr int] elki ' Rl7 ?i(Pi) $ i{ r 2}^xdridr 2 (15) 

with R j = ^(x+rj) and pj — (x — r j); j=l,2. Here, x is the coordinate of positron 
in Ps, and r^; j=l,2 are that of electrons in Ps and in H respectively in the incident 
channel w.r.t. the center of mass of the system. Functions r] and $ indicate the 
wave functions of Ps and H. Subscript ‘i‘ identifies the incident channel, whereas 
‘f represents the final channel. Accordingly k* and k f are the momenta of the 
projectile in the initial and final channels respectively. 


Vint = 



1 

I r l I 


1 1 

r "b i 

X - r 2 I | ri — r 2 


(16) 


V G - 

r nn 


/ int 



1 

I r 2 | 


1 

I x - ri | + 


1 

ri - r 2 | 


(17) 


Indices F and G indicate the direct and exchange channels respectively. 

The computation of Born-exchange matrix element is much more difficult than 
the direct Born-element. The Fourier transform 


the Bethe integral 


-A r 
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1 

J 

f e* pr 

27T 2 J 

p 2 + A 2 

Jk.r 

47T 

dr = — 

r 

k 2 


dp, 


and the properties of Dirac ^-function are mainly used for analytical evaluation. 
The numerical integrations are carried out using Gauss-Legendre quadratures. 


4 Results &; discussion 


In the present review article, we discuss the elastic scattering at low energies below 
the inelastic threshold of excitation in Ps-H and Ps-Li scattering. We use different 
target-inelastic and projectile- inelastic CCA schemes to study s-wave elastic phase 
shifts and s-wave elastic cross sections with a very fine mesh points to search for 
resonances- The incident energies are chosen from 0 to 5.1 eV for Ps-H system. 
Our results using the static-exchange approximation, target-inelastic 2-channel 
and 3-channel CCAs are presented in figures 1 & 2 and are giving no trace of 
resonances. But both the phase shifts and cross sections using projectile-inelastic 
2-channel and 3-channel CCAs are showing perfect resonances in singlet channel 
and are presented in figures 3 & 4. The singlet resonance using projectile-inelastic 
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Figure 7: Comparison of total cross sections for Ps-H scattering below inelastic 

thresholds using projectile-inelastic CCAs, static-exchange approximation and 22PslH 
coupled-pseudostate R-matix theory. 

3-channel CCA is shifted to lower energies and broadened in width than 2-channel 
CCA when it was very close to first excitation threshold. All the corresponding 
triplet results are also presented in the same figures 1 to 4. No triplet resonance 
is found except using the projectile-inelastic 3-channel CCA scheme. The triplet 
resonance is found near the energy 3.25 eV. The corresponding s-wave partial cross 
section is giving a sharp well. All resonances satisfy the Breit-Wigner formulation 
of cross section. This kind of perfect resonances were first time observed by Ray 
[1-2,68-72] using different projectile-inelastic CCA schemes. The singlet resonance 
was predicted earlier [25,35,46,66-67], but the triplet was new. 

We also reinvestigate the s-wave elastic phase shifts and cross sections below 
inelastic threshold for Ps-Li scattering using static-exchange approximation and 2- 
channel target-inelastic CCA in a similar fashion. Energy region is chosen from 0 to 
1.8 eV i.e. below inelastic threshold of first excitation. Our findings are presented 
in figures 5 & 6 and are very similar to the earlier results. The triplet channel is 
more sensitive to long-range forces due to dipole polarizabilty near to threshold. 

This long range dynamic effect [87-88] is introduced into our calculation through 
Li(2p) state. In a different way, it is the quantum mechanical effect of strong 
coupling between the Li(2s) and Li(2p) states. There is no trace of resonances 
in the singlet channel using both the CCA schemes. But the triplet channel is 
showing a rapid change in phase shifts ~ 7r/2 radian in both the static-exchange 
and 2-channel CCA at very close to zero energy and the corresponding partial wave 
cross sections are showing shaxp peaks. However these axe not resonances. The 
definition of below-threshold resonances admits a phase shift change of ir radian. 

It is not known to us what it indicates. 

In figure 7, we compare the total cross sections. These are again the total 
elastic cross sections since at the energy region below inelastic/excitation threshold 
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elastic channel is the only open channel. The projectile-inelastic 2-channel and 
3-channel CCA cross sections are presented with static-exchange approximation 
[47-48] and the 22PslH coupled-pseudostate R- matrix data of Blackwood et al [37]. 
The present results show good agreement. The well in the total cross section curve 
using projectile-inelastic 3-channel CCA supports the triplet resonance. The small 
peaks using projectile-inelastic 2-channel and 3-channel CCAs are supporting the 
singlet resonance. 

In addition, the present triplet phase shift data for Ps-H scattering using 
projectile-inelastic 3-channel CCA fit nicely with non-resonant part as 

£o = -1-4053 + 0.1295E - 0.04613E 2 


and provides the width T = 0.1517361^ and resonance position Er = 3.263061^. 

5 Conclusion 

A brief survey on Ps-atom scattering is made. Resonances below excitation thresh- 
old are discussed. Thorough reinvestigations are made in Ps-H system with static- 
exchange, target-inelastic and projectile-inelastic 2-channel and 3-channel CCAs 
using a very fine mesh points to search for resonances. No resonance is found with 
static-exchange approximation and target inelastic 2-channel and 3-channel CCAs 
in Ps-H system. A singlet resonance found in both projectile-inelastic 2-channel 
and 3-channel CCAs in Ps-H system which for the first time is detected using close- 
coupling approximation schemes, agree well with earlier predictions. A triplet res- 
onance with 3-channel projectile-inelastic CCA in Ps-H system is a new addition. 
It needs more investigations using an extended basis set and a systematic study 
to realize the physical cause of it. In Ps-Li system no resonance is found below 
inelastic threshold using both the static-exchange and target-inelastic 2-channel 
CCA. The comparison of total cross section curves below inelastic threshold in 
Ps-H system supports the existence of both the singlet and triplet resonances. 
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